Atom lithography with near-resonant hght masks: Quantum optimization analysis 
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We study the optimal focusing of two-level atoms with a near resonant standing wave light, using 
both classical and quantum treatments of the problem. Operation of the focusing setup is considered 
as a nonlinear spatial squeezing of atoms in the thin- and thick-lens regimes. It is found that the 
near-resonant standing wave focuses the atoms with a reduced background in comparison with far- 
detuned light fields. For some parameters, the quantum atomic distribution shows even better 
localization than the classical one. Spontaneous emission effects are included via the technique of 
quantum Monte Carlo wave function simulations. We investigate the extent to which non-adiabatic 
and spontaneous emission effects limit the achievable minimal size of the deposited structures. 

PACS number(s); 32.80.Lg, 03. 75. Be, 42.82. Cr, 03.65.Sq 



I. INTRODUCTION 

The ability to control the motion of atoms using laser 
fields has led to the realization of optical elements such 
as mirrors, lenses, beam splitters, etc, for atomic beams. 
One of the interesting applications is the laser focusing 
of atoms, which is useful to the technologically impor- 
tant problem of atom lithography. The principle of atom 
lithography is based on using a standing wave (SW) of 
light as a mask on atoms to concentrate the atomic flux 
periodically and create desired patterns at the nanome- 
ter scale (for recent reviews see, [1,2]). Since the first 
experimental demonstration of submicron atomic struc- 
tures [3], the subject has seen a considerable growth both 
theoretically [4-8] and experimentally [9-15]. In the di- 
rect deposition setup, periodic atomic lines of sodium 
[3,9], chromium [10,11], aluminum [12], ytterbium [13], 
and iron [14] atoms have been successfully fabricated. 
The technique has also been applied to two-dimensional 
pattern formation [15]. 

Most theoretical studies of atom lithography employ 
a particle optics approach to laser focusing of atoms 
[4,5,16]. The focal properties of the light have been ex- 
amined in terms of time dependent classical trajectories 
of atoms in the light induced potential. It has been shown 
that the atomic image at the focal plane exhibits a broad- 
ening due to severe aberrations caused by anharmonic- 
ity of the sinusoidal dipole potential [4,5]. As a result, 
all current lithography schemes suffer from a consider- 
able background in the deposited structures. A novel 
scheme to reduce the aberration problem was suggested 
recently in [16] by using optimized multilayer light masks. 
Quantum mechanical analysis of the focusing of atomic 
beams has been performed as well. Cohen et al. studied 
quantum mechanically the thin-lens regime of atom fo- 
cusing with both far detuned and exactly resonant stand- 
ing light waves [6] . Recently, atomic nanostructures have 
been experimentally realized with exactly resonant SW 
[17]. In the thick- lens regime, focusing atoms is generally 
achieved with a blue-detuned SW light whose detuning 
(A) is of the order of the Rabi frequency (ilg) of the 



atom-light interaction [4,5]. In this case, the influence 
of spontaneous emission on the focusing of atoms has 
been shown to be negligible by the quantum treatment 
[8]. However, a detailed study of the atom focusing with 
SW's that optimizes the lens performance as a function 
of the detuning of the light frequency has not appeared 
previously, to the best of our knowledge. 

In this paper, we study atom lithography with near- 
resonant light masks (|A|/f2o ^ 1) considering both blue- 
and red-detuned light. We give a comprehensive theoret- 
ical analysis for the classical and quantum treatments of 
the problem. The parameters for the best focusing of 
atoms are found as functions of the detuning (|A|) using 
the optimal squeezing method developed in Ref. [16]. To 
include momentum diffusion and spontaneous emission 
effects on atoms, we perform Monte Carlo wave func- 
tion simulations. High resolution deposition of chromium 
atoms is considered as an example, though the general 
conclusions drawn should apply well to other atoms. 

The paper is arranged as follows. In Sec. II, the ba- 
sic framework of the problem is defined and the focus- 
ing of atoms is studied classically under the influence 
of adiabatic light potentials. In this section, we exam- 
ine the optimal squeezing scheme of [16] when applied 
to the atomic-beam traversing a near-resonant SW light. 
In Sec. HI, the problem is treated in the quantum do- 
main to account for non-adiabatic and diffraction effects 
on atom focusing. The effects of spontaneous emission of 
atoms are considered in Sec. IV. Finally, in Sec. V, we 
summarize our main results. 



II. FOCUSING OF ATOMS BY ADIABATIC 
LIGHT POTENTIALS: CLASSICAL 
TREATMENT 

We begin our discussion with a model based on the 
interaction of a beam of two-level atoms with a near- 
resonant SW light. We take the direction of propaga- 
tion of the atomic beam through the light along the z 
direction. The SW (assumed to be formed along the x 
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direction) has a frequency of uji and A = tui ~ loq de- 
fines the detuning of the light frequency from the atomic 
transition frequency loq. The atom-light interaction is 
characterized by the Rabi frequency 



f2(a;, z) = r^o exp(— z^/cr^) cos(fca;) 



(1) 



Here, the term exp(— z^/tr^) accounts for the spatial vari- 
ation (a Gaussian beam profile with diameter a^) of the 
light intensity along the z direction. The cos{kx) term 
comes from the sinusoidal [I(x) oc cos^(fca;)] variation of 
the SW intensity along the x direction. The quantity 
ilo represents the peak Rabi frequency of the atom-light 
interaction and A = 27r/fc is the wavelength of the laser 
beams forming the SW. The velocity Vz of the atoms 
along the beam axis is sufficiently large, so that the 
atom's position along the z direction can be replaced by 
the time dependence z = v^t. Defining = cr^/vz-, the 
time dependent Rabi frequency is thus given by 



fl{x,t) — r^o exp(— t^/cTj ) cos(fca;) . 



(2) 



The behavior of atoms in the near-resonant light field 
can be best understood in the dressed state picture of the 
atom-light interaction [18]. The dressed states, which are 
the eigenstates of the interaction Hamiltonian, depend 
on time through the time-dependent Rabi coupling Eq. 
(2). If the Hamiltonian temporal variation is smooth, the 
atom prepared initially in one of the eigenstates of the 
Hamiltonian will follow the time-dependent eigenstate. 
The corresponding adiabatic condition is [19] 



A » J — 



(3) 



We assume that the atoms in the beam are initially in 
their ground state and that the adiabatic condition Eq. 
(3) is satisfied. In this case, the atoms can be described 
as pointlike particles moving in the potential 



U{x,t) 



h 



SgniA)-y/A-' + Q^x,t) 



(4) 



Here Sgn{A) = +1 (-1) for A > (A < 0). Note 
that this adiabatic potential should be contrasted with 
the light shift felt by the bare atomic states in the far- 
detuned limit [6]. The adiabatic potentials and the fo- 
cusing of atoms are shown schematically in Fig. 1. 

Many aspects of atom focusing by the standing light 
waves can be explained in the semiclassical picture of 
the atom's interaction with the light induced potential 
[4,5,16]. Therefore, we start considering the problem in 
the classical framework. We neglect spontaneous emis- 
sions from the atoms by assuming that the atom's in- 
teraction time with the light is much shorter compared 
to the lifetime of the excited atomic level. The classical 
trajectories of atoms in the adiabatic light potential Eq. 
(4) obey Newton's equation of motion 
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FIG. 1. (Color online) Schematic representation of laser fo- 
cusing of atoms by a SW light, (a) A coUimated atomic beam 
impinges on the near-resonant SW light whose intensity varies 
sinusoidally along the x direction, (b) The interaction with 
the hght induces the adiabatic potential [either plus or minus 
in Eq. (4)] for the external motion of atoms depending upon 
the sign of the detuning. For a positive (blue) detuning, the 
atoms localize near the minima of the light intensity, while for 
a negative (red) detuning they remain localized near the light 
intensity maxima. The dashed curves represent the quadratic 
approximation to the adiabatic potentials near their minima. 



where m is the atomic mass. The force induced by the 
SW focuses (localizes) the atoms near the intensity min- 
ima (maxima) for A > (A < 0). As a measure of the 
atomic localization, we use the localization factor [16] 



L{t) = 1 + Sgn{A){coa[2kx{t, xq)]) 
2 /•A/4 



A 



A/4 



dxo{l + Sgn{A)cos[2kx{t,xo)]} , (6) 



where x{t,xo) is the solution of the differential equation 
(5) satisfying the initial condition x —^ xq at t —^ — oo. 
The average in Eq. (6) is taken over the random ini- 
tial positions of the atoms, and the localization factor is 
measured as a function of time t counted from the mo- 
ment {t — 0) of passing the Gaussian center of the SW. 
The localization factor equals zero for an ideally localized 
atomic ensemble and is proportional to the mean-square 
variation of the x coordinate (modulo standing-wave pe- 
riod) in the case of a well- localized distribution (L ^ 1). 

The localization factor defined in Eq. (6) measures 
a nonlinear spatial focusing of atoms beyond the linear 
paraxial approximation. The plane of the best atomic 
localization (minimal atomic background) is determined 
by minimizing the localization factor [16]. This optimiza- 
tion can be done either in the thin- or the thick-lens limit 
of atom focusing. First we consider the thin-lens focusing 
[20] of atoms by the adiabatic potential Eq. (4), which 
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FIG. 2. (Color online) (a) Localization factor of the atomic 
distribution as a function of the dimensionless time t for 
at = 0.07 and A/Qq = -0.125 (solid curve), A/^o = 0.125 
(dashed curve), A/Qq = 5 (dotted curve). The global mini- 
mum values of L{t) are 0.17 (solid curve), 0.63 (dashed curve), 
and 0.42 (dotted curve), (b) Probability density (P) of the 
atomic distribution at the time t = tm of the best atomic 
localization. The parameters are same as those of (a) with 
tm = 0.82 (solid curve), tm ~ 0.52 (dashed curve), and 
tm = 7.38 (dotted curve). For the sake of comparison, the 
solid curve has been displaced by 0.25 units along the X axis. 
The times tm correspond to the global minima of the local- 
ization factor in (a). 



is valid in the Raman-Nath approximation to the atom- 
light interaction. In this case, the atomic displacement 
along the SW direction is negligible within the light and 
the focal point is well outside the region of the light fields 
[5]. Within thin-lens approximation, the SW introduces 
a spatially dependent sudden kick for atoms along the 
transverse {x axis) direction. For atoms with zero initial 
velocity along the x axis, the change in velocity can be 
calculated from Eq. (5) to be 



dvx = Sgn{A) 



4m 



dt- 



exp(-2t2/cr2) sin[2kx{t)] 
v/A2 + nlexp{~2t^/a^)cos^[kx{t)] 



(7) 



where the integration is done over the duration (±3crt) of 
the Gaussian profile. In the Raman-Nath approximation, 
the position x(t) of the atom in Eq. (7) during interac- 
tion with the light can be approximated to be its initial 
position xq. After passing through the light region, the 



0.7 



0.5 



0.1 





^^^^ (a) 







-5.0 -3.0 -1.0 1.0 3.0 5.0 

A/n„ 

Thick Lens - Classical 




FIG. 3. (Color online) Minimal localization factor of the 
atomic distribution as a function of the detuning A/Qq for 
(a) at = 0.07, (b) at = 4. 

atom moves as a free particle with a time dependent 
transverse position 

x{t,xo) = xo + Svxt . (8) 

Using Eqs. (6) and (8), we calculate the localization 
factor as a function of dimensionless time for both blue 
(A > 0) and red (A < 0) detuning conditions. The 
results are shown in Fig. 2 with plots of the atomic dis- 
tribution [21] at the time of minimal localization factor. 
These results are also compared with the best atomic 
localization that can be achieved with the usual far- 
detuned (A/r^o > 1) SW hght. We use dimension- 
less variables in which position is measured in units of 
the optical wavelength (A), frequency in units of the re- 
coil frequency {oJrec = hk'^/2m), and time in units of 
l/[ujrec^o'^t]- It is seen from the graph that the near- 
resonant light with red detuning from the atomic tran- 
sition focuses the atoms better compared with the far- 
detuned light. The atomic localization improves by de- 
creasing the magnitude of the detuning in the case of 
atom focusing by the red-detuned light. This is shown in 
Fig. 3(a), where we plot the best localization factor as 
a function of the detuning parameter (A/ilo). Note that 
the localization factor saturates to the well known value 
of L « 0.42 in the far-detuned (|A|/17o > 1) case [16]. It 
tends to the asymptotic value of L « 0.1 for A ^ 0^. We 
note, however, that nonadiabatic quantum effects should 
be taken into account when the adiabatic condition, Eq. 
(3) fails (see the next section). 
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FIG. 4. (Color online) (a) Localization factor of the atomic 
distribution as a function of the dimensionless time t for 
at = 4, and A/fJo = -0.125 (sohd curve), A/ilo = 0.125 
(dashed curve), A/Qg — 1 (dotted curve). The global mini- 
mum values of L{t) are 0.1 (solid curve), 0.36 (dashed curve), 
and 0.31 (dotted curve), (b) Probability density (P) of the 
atomic distribution at the time t — tm of the best atomic 
localization. The parameters are same as those of (a) with 
tm = 0.2 (solid curve), tm — 0.58 (dashed curve), and tm ~ 2 
(dotted curve). For the sake of comparison, the solid curve 
has been displaced by 0.25 units along the X axis. The times 
tm correspond to the global minima of the localization factor 
in (a). 



In the opposite limit of the thin lens, the atoms are fo- 
cused within the region of the light fields [5] . The atom 
focusing in this limit (known as the thick-lens focusing) is 
very similar to the operation of the graded index lens in 
traditional optics. Neglecting the time dependence in the 
Rabi frequency, we find the focal length in the paraxial 
approximation as 



2^0 V 



[A > 0, Thick lens] , 



f=^\ [A < 0, Thick lens] . (9) 

From the above expression, it is seen that the focal length 
is greater for red detuning. To go beyond the paraxial 
approximation, the solution x(t) of the atomic motion 
needs to be obtained directly by the numerical integra- 
tion of Eq. (5). Using Eqs. (5) and (6), we obtain the 
localization factor of the atomic distribution as shown in 
Figs. 3(b) and 4. 



On comparing the graphs from Figs. (2) to (4), we 
see a qualitative similarity between the thin- and thick- 
lens focusing of atoms. In particular, focusing atoms by 
the red-detuned light gives rise to a much reduced back- 
ground of deposited atoms in comparison to the case of 
atom focusing by the blue-detuned light. In the paraxial 
approximation, this can be explained by expanding the 
adiabatic potentials Eq. (4) upto the quadratic terms 
(parabolic fitting) near the potential minima [22]. As 
shown in Fig. 1, the spatial range, in which the quadratic 
approximation is valid, is wider for the case of negative 
detuning (A), which results in the reduced aberrations. 
We note that the above analysis is valid only in the range 
of parameters for which the adiabatic condition Eq. (3) 
is satisfied. For smaller detunings (|A|/r2o ^ 1), non- 
adiabatic and quantum effects may dominate which will 
be examined in the next section. 



III. OPTIMAL ATOMIC SQUEEZING - EFFECTS 
OF THE WAVE NATURE OF ATOMS 



In the classical treatment of atom focusing discussed 
so far, the internal structure of the atom was completely 
ignored and the atomic motion in single adiabatic po- 
tential was studied in the time domain. This procedure 
is valid even if the atomic motion is treated quantum 
mechanically, provided the adiabatic condition (3) is sat- 
isfied. In the quantum treatment, the atomic center-of- 
mass (cm.) wave function tp^x^t) evolves in time ac- 
cording to a Schrodinger equation in which the potential 
U{x,t) is given by Eq. (4): 



d 



^ + U{x,t) 
2m 



il){x,t) , 



(10) 



where denotes the cm. momentum operator of the 
atom along the SW {x axis) direction. 

While the time evolution, Eq. (10) is useful to study 
the quantum effects on atom focusing, the situation be- 
comes more complicated if the light detuning is relatively 
small (|A| <C r^o)- In this case, non-adiabatic effects 
arise from transitions between the dressed atomic states. 
Therefore, in order to cover uniformly a wide range of 
detunings, we consider the evolution of the atomic wave 
function directly in the bare-states basis, first neglecting 
spontaneous emission from the atoms. The Hamiltonian 
for a two-level atom with excited (|e)) and ground (|.g)) 
states interacting with the SW light is given by 



Hit) 



2m 2 
^hU{x,t) 



i\e){e\ 
{\e){9\ 



\9){9\) 
\9){e\) 



(11) 



where f2(x,t) is given by Eq. (2). 

f,,„„<-;^„ ^]-^g two-level atom may be ex- 
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FIG. 5. (Color online) Localization factor of the atomic 
distribution as a function of the dimensionless time t for 
the parameters at = 0.0006, Qo = 1-92 x 10^, and (a) 
A/ilo = -0.125, (b) A/ilo = 0.125 (solid curve), A/^o = 5 
(long-dashed curve). The global minimum values of L{t) are 
(a) 0.15, and (b) 0.2 (solid curve), 0.42 (long-dashed curve). 



(12) 



Here il}(,,g{x,t) correspond to the cm. wave functions of 
the atom in its excited and ground states. We consider a 
spatially uniform beam of ground state atoms having ini- 
tially zero momentum along the SW direction. The ini- 
tial wave function of the atom normalized over the region 
of the SW period [21] is then given by il^g{x, to) = ^2/A, 
where the initial time to — > — oo. Since the atomic distri- 
bution is expected to be periodic (in space) after inter- 
action with the light field, the wave functions at time t 
can be Fourier expanded as 



n— — oo 
n— oo 

i^g{x,t)^ J2 Cme''"'^ ■ (13) 
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FIG. 6. (Color online) Probability density (P = |*(j:,t)p) 
of the atomic distribution at the time t = tm of the best 
atomic localization. The parameters are at — 0.0006, 
= 1.92 X 10^ and A/!^o = -0.125, t„ = 7 x 10"^ (solid 
curve), A/r2o = 0.125, = 0.778 (dashed curve), A/Qq = 5, 
tm = 0.72 (long-dashed curve). For the sake of comparison, 
the solid curve has been displaced by 0.25 units along the X 
axis. The times tm correspond to the global minima of the 
localization factor in Fig. 5. 



'jC:Xt) = l{2n+irL0r.,-^}CUt) 



n{t) 



(14) 



n{t) 



(Q(i) + Q_i(i)) 



with n{t) = f^oexp(-i2/cr2). 

We now proceed to the calculation of the quantum lo- 
calization factor defined as 



L{t) = 1 + S',gn( A) (cos(2fca;)) 
rX/4 



l + Sgn{A) / dx\'l!{x,t)\ cos{2kx) 



(15) 



The atomic density 1^(0:, i) I —\ipe{x,t)\ +\%lig{x,t)\ is 
found using Eqs. (13). For the initial beam of ground 
state atoms (Cq(— 00) = a/2/A), the localization factor 
can be obtained by solving numerically Eqs. (14): 



L{t) = 1 + S'gn(A) 



A 



n— — 00 



;Re [Cn{t)C'n+i{tr + Cl{t)Ci^,{tn (16) 



The Fourier coefficients C^(t) and C^{t), defined above, 
represent the probability amplitudes for finding the atom 
in the excited and ground states with momentum (2n + 
l)hk and 2nhk, respectively. Using Eqs. (11)-(13), the 
Schrodinger equation ihd'i//dt = then leads to cou- 
pled equations for the Fourier amplitudes as 



As discussed in the previous section, the best spatial 
squeezing (localization) of atoms occurs at the time 
t = tm of the global minimum of the localization fac- 
tor. To minimize the localization factor Eq. (16), we 
first rescale the time variable t in terms of the recoil time 
tree = In the classical analysis of atom focusing, 



5 



Thin Lens - Quantum 



0.5 



0.4 



L(tJ 0.3 



0.2 - 



0.1 



(a) : 







-5.0 -3.0 -1.0 1.0 3.0 5.0 
Thick Lens - Quantum 




FIG. 7. (Color online) Minimal localization factor of the 
atomic distribution as a function of the detuning A/f2o for 
(a) at = 0.0006, (b) at = 0.01. 

we used a time scale that depends on the strength (fio) of 
the atom-light interaction. However, a quantized space- 
periodic motion of ground state atoms in free space 
[23] repeats itself (in time) after each revival period of 
tfi = 7rtrec/2. It is therefore convenient to use the recoil 
time tree for scaling the variables t and at- 

In Figs. 5 and 6, we display the localization factor and 
the atomic spatial density for the thin-lens focusing of 
atoms. On comparing the quantum mechanical localiza- 
tion factor with its classical counterpart (see Figs. 2 and 
5), we see that they have the same structure for short 
times and the same minimal values for both blue and 
red detuning situations. However, the localization fac- 
tor differs considerably in the long time evolution of the 
atomic distribution. It reaches a new minimum value of 
L — 0.2 at about half of the revival period for the blue 
detuned {A/Hq = 0.125) light. Due to these distinctly 
quantum features in the atom localization, the quantum 
atomic distribution becomes narrower than the classical 
distribution (compare Figs. 2(b) and 6) in the atom fo- 
cusing by the blue detuned light. Next, we show the best 
quantum localization of atoms that can be achieved by 
varying the detuning in Fig. 7(a). Again comparing it 
with the classical result (see Fig. 3(a)), it is seen that 
they are identical for the red detuning (A < 0) case. New 
quantum features exist only in the focusing of atoms by 
the blue detuned light (A > 0). For the thick-lens focus- 
ing of atoms, the results are shown in Figs. 7(b) and (8). 
In this case, the quantum and classical results are quite 



L(t) 



1 

0.8 

0.6 
0.4 
0.2 



: (a) 

'- 


\\ 

' / 


• \ 1 ' • r 

- V/ ' : 
\ ■ 




W 

\ r 




1 1 1 


1 





f.05 



0.05 




0.35 



FIG. 8. (Color online) (a) Localization factor of the atomic 
distribution as a function of the dimensionless time t for the 
parameters at = 0.01, r2o = 4 x lO", and A/fio = -0.125 
(solid curve), A/Slo ~ 0.125 (dashed curve), A/Qo ~ 1 
(long-dashed curve). The global minimum values of L{t) are 
0.1 (solid curve), 0.37 (dashed curve), and 0.31 (long-dashed 
curve), (b) Probability density (P = |'I'(a;, t) |^) of the atomic 
distribution at the time t = tm of the best atomic localization. 
The parameters are same as those of (a) with tm = (solid 
curve), tm ~ 1.5 x 10~^ (dashed curve), and tm ~ 5.3 x 10"'^ 
(long-dashed curve). For the sake of comparison, we have 
displaced the solid curve by 0.25 units along the X axis and 
the dashed curve by 20 units along the P axis. The times tm 
correspond to the global minima of the localization factor in 
(a). 



similar in the adiabatic regime of atom focusing and there 
are no distinct quantum features in the localization of 
atoms. The wave effects in propagation have a limited 
manifestation in this case due to the relatively short focal 
length. 

When the light frequency is tuned very close to the 
atomic resonance (|A| <C ^o), the adiabatic potentials, 
Eq. (4) experience sharp spatial variations along the x 
coordinate in the regions of avoided crossings near the 
nodes of the SW. As a result, non-adiabatic effects de- 
grade greatly the focusing of atoms. In this case, the 
atoms do not follow single adiabatic potential (either plus 
or minus in Eq. (4)) during interaction with the SW but 
rather make random transitions between them in the re- 
gions of quasi-crossing. As a result, they tend to focus 
near both the minima and maxima of the light intensity 
with no well defined localization region as shown in Fig. 
9. The localization factor defined in Eq. (15) is not 
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FIG. 9. Probability density (P = |*(a:,t)| ) of the atomic 
distribution for the parameters at = 0.01, f2o = 2 x 10*, t = 0, 
and A/fio = 10"^ 



appropriate to characterize the atomic distribution for 
this case. 



IV. MONTE CARLO SIMULATIONS 

We have so far simplified the classical and quantum 
treatments of the problem by assuming that the atoms 
do not change their internal states by spontaneous emis- 
sions diu'ing interaction with the light. This is strictly 
valid only if Tat ^ 1, where T is the decay rate of the ex- 
cited atomic level. However, for chromium atoms which 
we consider, the quantity Fctj becomes of the order of 0.1 
and 2 in the thin {at = 0.0006) and thick {at = 0.01) lens 
regimes, respectively. Therefore, it is necessary to study 
the extent to which spontaneous emission degrades the 
focusing performance of the adiabatic light potentials. 

The effects of spontaneous emission on the focusing of 
atoms are two- fold. First, it may modiiy the trajectories 
of atoms by imparting momentum kicks in random di- 
rections. Secondly, the dressed internal state of the atom 
changes due to spontaneous emissions resulting in ran- 
dom jumps between the plus and minus adiabatic poten- 
tials, Eq. (4) which the atom experiences. To model all 
these features, we employ the quantum Monte Carlo wave 
function (MCWP) simulations [24,25] for the description 
of atoms emitting spontaneous photons and moving in 
a SW light field. In this approach, the atomic density 
P is calculated by forming N realizations of quantum 
trajectories ipl^g{x,t) and then averaging over them : 

1 ^ 

P^\^{x,t)f = -^M{x,tf + \r,{x,t)f] . (17) 

s=l 

Each quantum trajectory ipl_g{x,t) consists of a set of 
deterministic Hamiltonian evolution periods interrupted 
by quantum collapses. The Hamiltonian evolution is gov- 
erned by a non-Hermitian Hamiltonian i?e// = H{t) — 
i?i(r/2)|e)(e|, where the term H{t) describes the coher- 
ent dynamics of the atom-light interaction as given in Eq. 
(11) and the imaginary term (containing F) accounts for 



the spontaneous decay of the excited atomic level. The 
collapse of the atomic wave function is given by the ac- 
tion of the operator 



= [TN{k'f^eM-ik'x) \g){e\ 



where 



N{k') = (3/8fc) 



-4 



(18) 



(19) 



is the normalized probability density for the distribution 
of the spontaneously emitted photons with momentum 
component Tik' along the SW direction {x axis). 

When there is no spontaneous emission event, the in- 
ternal states of the atom are coupled only by stimulated 
processes. The atomic wave function at time t is given 
by a Fourier series [cf. Eq. (13)] 

n—oo 

n= — oo 

n—oo 

ijg{x,t)= C9(i,ti) e*[f°+2"'^'=l^/'^ . (20) 



where po is the momentum of the atom along the SW 
direction at the initial time ti. The Fourier coefficients 
defined above depend on the initial time ti and evolve 
with time t (until a spontaneous emission takes place) as 
governed by the non-Hermitian Hamiltonian Hgf f : 



'Jt^^^' - 2^^ 



A 



cm 



(21) 



If a spontaneous emission from the atom takes place at 
the time t, the momentum hk' of the spontaneously emit- 
ted photon is chosen randomly according to the proba- 
bility law N{k') [Eq. (19)] and the collapse of the atomic 
wave function to the ground state is carried out with the 
operator Eq. (18) as follows 



cm = , 

Po ^ Po + ^k — hk' . 



(22) 
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L(t) 0. 



0.015 




FIG. 10. (Color online) (a) Localization factor of the 
atomic distribution as a function of the dimensionless time 
t around its global minimum for the parameters at = 0.0006, 
no = 1.92 X 10^ A/Qo = -0.125, and T = 238 (solid curve), 
r = (dashed curve). The minimum values of L{t) are 0.175 
(solid curve) and 0.15 (dashed curve), (b) Probability density 
(P = |^'(2;,t)|^) of the atomic distribution at the time t = tm 
of the best atomic localization. The parameters are same as 
those of (a) with tm = 7 x 10~^. The solid curve is almost 
indistinguishable from the dashed curve on the scale shown. 
For clarity, the close-up to the right of origin {x = 0) is shown 
in the inset. 



In the MCWF simulations [25], the random moment 
t at which the spontaneous emission takes place is cho- 
sen when the decaying total norm of the atomic wave 
function reaches the value of 1 — e, where e G [0, ...1] is 
a random number uniformly distributed between and 
1. The moment of emission is determined by solving the 
equation 



n— oo 

l-(A/2) ^ [\C:{t,U)f + \C?,{t,UW 



e . 



We take the initial time ti to be —hut (instead of 
—oo) at which the Rabi frequency ^{x^€) [Eq. (2)] 
drops significantly compared to its peak value. We as- 
sume, as before, the initial condition in the form of 
an atomic plane wave in the ground state with zero 
transverse momentum (po — 0). The norm |^(i)|^ — 



0.04 




FIG. 11. (Color online) (a) Localization factor of the 
atomic distribution as a function of the dimensionless time 
t around its global minimum for the parameters at = 0.01, 
fio = 4 X 10*, A/^o = -0.125, and F = 238 (solid curve), 
F = (dashed curve). The minimum values of Lit) are 0.18 
(solid curve) and 0.1 (dashed curve), (b) Probability density 
(P = |^'(x,t)|^) of the atomic distribution at the time t — tm 
of the best atomic localization. The parameters are same as 
those of (a) with t,„ = 0. For clarity, the close-up to the right 
of origin {x = 0) is shown in the inset. 



(A/2)E„ mt,u)r + mt,uw 



of the wave function 



is then obtained by solving Eqs. (21) until the time 
t ~ thic at which Eq. (23) is fulfilled. At this moment, 
we collapse the atomic wave function to the ground state 
and add a recoil momentum to the atom [see Eqs. (18) 
and (22)]. After the collapse, the random number e is 
renewed and the process gets repeated by solving again 
Eq. (23) starting from the new initial time tt = Imc and 
the new values for the initial Fourier coefficients. This 
procedure gets continued until the final time t at which 
the calculation of the deposited atomic density should 
be performed. This results in single quantum trajectory 
(23) -01 obtained by normalizing the wave function Eqs. 

(20) at the final time t. An ensemble average of many 
such trajectories [Eq. (17)] is statistically equivalent to 
the solution of the density matrix equations. 

In the limit, when the detuning of the light frequency is 
relatively large (|A| ^ 17o), the atoms primarily evolve 
in their ground state during interaction with the light. 
In this case, the degradation of the atom focusing due 
to spontaneous emission is negligible as reported in ear- 
lier studies [6,16]. However, if a near-resonant light is 
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used to focus the atomic beam, the spontaneous emis- 
sions may broaden significantly the atomic distribution 
because of the random atomic recoils and fluctuations of 
the focusing potentials discussed above. These effects are 
especially important in the atom focusing by red-detuned 
light fields, since, in this case, the focusing occurs near 
the light intensity maxima resulting in substantial pop- 
ulation of excited atoms. In order to measure this quan- 
titatively, we again use the localization factor Eq. (15) 
with the atomic density |\E'(x,t)| given by Eq. (17). We 
find that an average oi N = 5000 quantum trajectories is 
sufficient to give a statistical error below 2% of the mean 
value in Eq. (17). The results are shown with (F = 238) 
and without (F = 0) spontaneous decay of the atoms for 
the thin- lens (Fig. 10) and thick- lens (Fig. 11) focus- 
ing regimes. It is seen from the graphs that spontaneous 
emission does not shift the optimal time t = tm for the 
global minimum of the localization factor. But, it re- 
duces the atomic density and increases the background 
significantly in the thick-lens focusing of atoms. 

V. SUMMARY 

In this paper, we studied the spatial focusing of an 
atomic beam by a near-resonant SW light in the con- 
text of atom lithography. The problem was treated both 
classically and quantum-mechanically using the optimal 
squeezing approach suggested recently [16] for the effi- 
cient focusing of atomic beams. Both thin- and thick- 
lens regimes were considered with blue- as well as red- 
detuned light fields. In the thick-lens regime, we have 
shown that a red-detuned light focuses the atoms bet- 
ter in comparison with a blue-detuned one and the op- 
timal quantum localization of atoms follows the classical 
scenario. In the thin-lens regime, we have shown that 
quantum-mechanical effects give rise to new localization 
mechanisms which are more effective than the (Jassical 
ones. Finally, we considered the role of non-adiabatic 
transitions between the dressed atomic states and the ef- 
fects of spontaneous emission of atoms on the quality of 
the deposition profile. 
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